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A relativistic theory of modified gravity has been recently proposed by Bekenstein. The tensor 
field in Einstein's theory of gravity is replaced by a scalar, a vector, and a tensor field which interact 
in such a way to give Modified Newtonian Dynamics (MOND) in the weak-field non-relativistic limit. 
We study the evolution of the universe in such a theory, identifying its key properties and comparing 
it with the standard cosmology obtained in Einstein gravity. The evolution of the scalar field is akin 
to that of tracker quintessence fields. We expand the theory to linear order to find the evolution 
of perturbations on large scales. The impact on galaxy distributions and the cosmic microwave 
background is calculated in detail. We show that it may be possible to reproduce observations of 
the cosmic microwave background and galaxy distributions with Bekenstein's theory of MOND. 

PACS numbers: 98.90.Cq, 98.65.Dx, 98.70.Vc 
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The current model of the Universe is based on a few 
simple assumptions and can explain a multitude of ob- 
servation. Yet, to be able to explain the structure of 
galaxies and clusters of galaxies, it is essential to postu- 
late the existence of some invisible substance, called dark 
matter. Although there are reports of tentative discover- 
ies of dark matter p| , there is no proven theory or direct 
observation of a dark matter particle as yet. A less ex- 
plored route is that our current theory of gravity might 
be incomplete. Given that the existence of dark matter 
is inferred from its gravitational effects on the dynamics 
of astrophysical bodies, it may be that the theory used to 
link the dynamics to the mass is incorrect. Milgrom has 
proposed a modification of Newtonian dynamics, known 
as MOND 0. In MOND, Newton's second law in a grav- 
itational field is modified to fi(\d\/ao)a = — V$ where a 
is the acceleration, $ is the Newtonian potential and fi(x) 
is a function with a scale x ~ 1. For \a\ > ao we have 
that n(x) = 1. With \a\ < ao we have that fi(x) ~ x. 
Clearly for small accelerations, Newtonian theory is no 
longer valid. Milgrom's theory has been extremely suc- 
cessful in explaining a number of observational properties 
of galaxies [3j. It has suffered from a fatal flaw in that 
it is not generally covariant and hence cannot be studied 
in a general setting. 

Bekenstein has recently solved this problem Q . Build- 
ing on a series of developments he has proposed a 
generally covariant theory which has in the nonrelativis- 
tic, weak- field limit, Milgrom's modified theory. Beken- 
stein's theory has two metrics. One of the metrics, 
has its dynamics governed by the Einstcin-Hilbert action, 
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where G is Newton's constant and R is the scalar cur- 
vature of g^ v . We shall call the frame of this metric 



the "Einstein Frame" (EF). The second metric, g^ is 
minimally coupled to all the matter fields in the Uni- 
verse. We shall call the frame of this metric the "Matter 
Frame" (MF). All geodesies are calculated in terms of 
this second metric. The two metrics are related through 
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Two fields are re- 



quired to connect the two metrics. The scalar field, 
dynamics given by the action 
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where /x is a nondynamical field and V is a free function 
which can be chosen to give the correct nonrelativistic 
MOND limit and depends on two free parameters, £b and 
Ho (related to k in 4] as fio = 8tt/k). The unit timelike 
vector field, An has dynamics given by the action 
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where F^ v — A^ tlJ — A v>fi , indices are raised with g and 
where K is the third parameter in this theory. The La- 
grange multiplier A is completely fixed by variation of the 
action. 

We wish to study the evolution of a homogeneous and 
isotropic universe in such a theory. Observers are defined 
in the MF where the line element is ds 2 — a 2 (~drj 2 -\-dr 2 ). 
The scale factor is related to the metric in the EF, ds 2 = 
b 2r_ e -4<f, dr] 2 + dr 2) though a = be - *. The modified 
Friedmann equation becomes 
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where fi can be found by inverting <f> 2 — \a 2 e 2 *^j and 
the energy density p does not include <f>. Homogeneity 
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and isotropy and the constraint in the action imply that 
An is fixed as A„ = ae _ ^(l, 0, 0, 0). The background 
dynamics is complete with an equation for <fi: 



-a 2 e~ 



2<j>yl 



1 

U 



2(p 



V" ~t> ( 



+47rGa 2 e- 4 ^( j o + 3P) 



where U = /x + 2V' /V" and, once again, P does not 
include the pressure from the scalar field 4>. 

A choice of V will pick out a given theory. As a first 
guess, Bekenstein has proposed 
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where ft = p/pa- This potential will lead to the pre- 
scription proposed by Milgrom in the nonrelativistic re- 
gion. The evolution of this coupled set of equations is 
essentially insensitive to £b and independent of K but is 
well determined in terms of po- If we define the physical 
Hubble parameter H, we can rewrite the modified Fried- 
mann equations in the form 2>H 2 = d>irG e ff{p+p ( f ) ) where 
the effective Newton's constant has the form G e // = 

Ge~ 4, 7(l + rftrs;) 2 and the ener gy density in <p is 
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This system exhibits the tracking behavior witnessed in 
some scalar field theories of quintessence 0. In particu- 
lar one finds that for a wide range of initial conditions, 
(j) evolves to a slowly varying function of time and the 
relative energy density in (f> reaches an attractor solu- 
tion of the form Jl^ = 3/(2/io) in the radiation era and 
fi^ = l/(6/io) in the matter and A eras. 

In Fig. 1 we solve the equations numerically for a low 
value of po to illustrate the tracking behavior of 4>. As 
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FIG. 1: The relative energy densities in <f> (thick solid line), ra- 
diation (dotted line), matter (dashed line) and A (dot-dashed 
line) for /io = 5 as a function of the scale factor (a is in arbi- 
trary units). Note that the energy density in the scalar field 
tracks the dominant form of energy at each instance in time. 



density in at nucleosynthesis 0. The abundance of 
light elements is extremely sensitive to the expansion rate 
at 1 MeV and can, for example, be used to constrain the 
number of relativistic degrees of freedom at that time. 
Recent measurements of the 4 He mass fraction and the 
deuterium abundance leads to a bound at energies of 1 
MeV of flcf, < 0.045 at the 95% confidence level. This 
leads to po > 33 and hence the scalar field will make 
up less than 0.5% of the total energy density during the 
matter and A dominated eras. 

We now turn to evolution of linear perturbations on 
this background. This will allow us to link this theory 
with observations of galaxy clustering on large scales as 
well as with the anisotropies in the Cosmic Microwave 
Background (CMB). The main problem that a MOND 
theory containing only baryons has to confront is the 
damping of perturbations during the recombination era. 
Indeed in a pure baryonic universe evolving under Ein- 
stein gravity, the weak coupling of baryons and photons 
during the recombination era will lead to Silk damping, 
the collisional propagation of radiation from overdense to 
underdense regions 8]. In the standard adiabatic model 
with just baryons, the matter power spectrum is severely 
suppressed on galactic scales. If MOND is to succeed, it 
must overcome the Silk damping on these scales. 

There are three fields to perturb around the back- 
ground (we identify the homogeneous part of a quantity 
X by X). The MF metric perturbations are given in 
the Conformal Newtonian gauge by goo = — a 2 (l + 2 1 i>) 
and gij — a 2 (l — 2Q)5ij. The scalar field perturba- 
tion is given by 4> — 4> + <P- The vector field pertur- 
bation is defined by — ae^^(A fl + ct M ) where the 
scalar components, a and E, of a M = — ip, a) are 
given by V 2 a = V • a and V 2 E = V • E where we 
use the field strength tensor of A^, F M „ to define the 
"electric field" through Ei — A^F itl . Similarly, in EF, 
we have 300 = -b 2 e~ 4 ^(l + 2*), g Qi = -b 2 (.i and 
gij = b 2 (l - 2$)(S y - which give 9 = ^ — tp : $ = <f> — <p 
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The evolution equations for the matter fluid remain 
unaltered if expressed in terms of the MF variables. That 
is if we expand densities as p = p(l + S) and use the 
standard definition for momentum of the fluid, V 2 6> = V- 
v, the evolution equations remain the same as in Einstein 
gravity. Two new sets of evolution equations must be 
introduced. For the scalar field perturbations, we have 
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and for the vector field we have 
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with other tracking systems we can constrain the energy 
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where w = P/p and c 2 = 5P/5p. The perturbed Einstein 
equations allow us to identify the gravitational potentials 
through: 
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where a is the total shear from matter fluids. 

We have modified CMBEASY, a publicly available 
numerical Einstein-Boltzmann solver to incorporate the 
modified background and perturbation equations 
The evolution equations have been implemented in both 
the conformal Newtonian and synchronous gauge to 
check for consistency [lOj- We have restricted ourselves 
to a flat Universe with a cosmological constant but con- 
sidered the possibility of massive neutrinos. 
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FIG. 2: The effect of the MOND parameters on the power of 
spectrum of the CMB. Top panel: p = 200, £ B = lOOMpc 
and K = 1 (solid), 0.1 (dotted) and 0.08 (dashed); Middle 
panel: p = 200, K = 0.1 and l B = lOOOMpc (solid), lOOMpc 
(dotted) and lOMpc (dashed); Bottom panel: K = 0.1, £b = 
lOOMpc and p = 1000 (solid), 200 (dotted) and 150 (dashed). 

The parameters po, K and Ib may introduce major 
modifications in the morphology of perturbations. A low 
Po, low Ib and low K will lead to a change in the growth 
rate. As we can see in Fig. 2, the effect is to introduce 
an integrated Sachs Wolfe term which can be quite sig- 
nificant. For example, for sufficiently small £b, the struc- 
ture of the angular power spectrum of the CMB can be 



completely modified with an excess of large scale power 
overwhelming structure on the smallest scales. We can 
see the effect of modifying £b in the lower panel of Fig. 2. 
Clearly, the CMB can place quite stringent constraints on 
the values of these parameters. A further possible effect, 
that we have not included here, is the effect of secondary 
anisotropics such as the Sunyaev-Zel'dovich or Ostriker- 
Vishniac effects, which may leave a different signature 
than the standard cold dark matter ACDM model. 

Since the baryon content is set by the abundance of 
light elements, we must compensate with a high value 
of the cosmological constant, i.e. with fl\ ~ 0.95. An 
obvious consequence of this is that the angular-distance 
relation will be modified as compared to the standard 
adiabatic ACDM universe Indeed the position of 

the peaks in the angular power spectrum of the CMB 
will be shifted to higher Is which would lead to a se- 
vere mismatch with the current available data from the 
Wilkinson Microwave Anisotropy Probe and other exper- 
iments. A natural solution to this is to include a small 
component of massive neutrinos, £l v ~ 0.15. As we can 
see in the top panel of Fig. 4, with this modification we 
can reproduce the temperature anisotropy data. 
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FIG. 3: The effect of the MOND parameters on the power 
of spectrum of the baryonic density fluctuations. Top panel: 
Ho = 200, £ B = lOOMpc and K = 1 (solid), 0.1 (dotted) 
and 0.08 (dashed); Middle panel: p = 200, K = 0.1 and 
l B = lOOOMpc (solid), lOOMpc (dotted) and lOMpc (dashed); 
Bottom panel: K = 0.1, £ B = lOOMpc and ^ = 1000 (solid), 
200 (dotted) and 150 (dashed). 

The main question we have raised is whether MOND 
dynamics can inhibit the damping of small scale pertur- 
bations in the coupled baryon-photon fluid during recom- 
bination. Recall that in the adiabatic CDM model, per- 
turbations in the dark matter, 5c, are undamped dur- 
ing recombination. The Newtonian potential, which is 
roughly given by fc 2 $ ~ AttG(pbSb + Pc&c) will not be 
erased if pc is sufficiently large, even though 5b — > 
through recombination 12]. In the MOND universe we 
find an analogous effect; we now have A: 2 $ ~ 4itGpb{5b — 
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2tp). The perturbation in the scalar field will support the 
perturbations through recombination yet still allow the 
damping of anisotropies in the photon fluid. Unlike the 
case of dark matter however, the coupling between the 
scalar field and the metric is such that p§ does not play 
a role in the magnitude of the effect. Even for minute 
values of fl^ we can still have a non-negligible effect. As 
we can see in Fig. 3, the net result is that decreasing 
/io, Ib or K will boost small scale power in such a way 
as to overcome the damping of perturbations. This is an 
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FIG. 4: The angular power spectrum of the CMB (top panel) 
and the power spectrum of the baryon density (bottom panel) 
for a MOND universe (with a ~ 4.2 x 10~ 8 cm/s 2 ) with fi A = 
0.78 and Q v = 0.17 and fi s = 0.05 (solid line), for a MOND 
universe Qa = 0.95 and fi_g = 0.05 (dashed line) and for the 
A-CDM model (dotted line). A collection of data points from 
CMB experiments and Sloan are overplotted. 



intriguing effect that goes in tandem with what we saw in 
the CMB. While decreasing Ib (and a sufficiently small 
K and /io) will contaminate the large scale power in the 



angular power spectrum of the CMB, it can also play a 
role in counteracting Silk damping of density perturba- 
tions. 

Given these two effects on the dynamics of large scale 
structure, is it possible to construct a MOND universe 
which can reproduce current observations of the CMB 
and galaxy surveys? There is clearly a competition be- 
tween overproducing large scale power in the CMB but 
also overcoming damping on small scale. In Fi g. 4 we 
present two MOND universes compared to data |l3t Il4| . 
As mentioned above, a universe with a very large contri- 
bution of A will not fit the current CMB data. By having 
the three neutrinos with a mass of m v ~ 2 eV each we 
are able to resolve this mismatch. With an appropriate 
choice of K, (j,q and Ib it is possible to reproduce the 
power spectrum of galaxies as inferred from the Sloan 
Digital Sky Survey [14J . The possibility of using massive 
neutrinos to resolve some of the problems with clusters 
in a MOND universe has been mooted in [l5j . 

We have focused on one very specific model proposed 
by Bekenstein with a somewhat artificial potential for 
the new degrees of freedom. This phenomenological ap- 
proach needs a firmer theoretical underpinning which 
might come from the various approaches which are being 
taken in the context of brane worlds, M-theory and a rich 
array of theories of modified gravity. However, Beken- 
stein's theory can play an important role in opening up 
an altogether different approach to the dark matter prob- 
lem. It serves as a proof of concept which will clearly 
lead to a new, very different view of the role played by 
the gravitational field in cosmology. 
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